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Abstract: We study on the lattice the topology of SU(2) and SU(3) Yang-Mills theories at zero tem- 
perature and of QCD at temperatures around the phase transition. To smooth out dislocations and 
the UV noise we cool the configurations with an action which has scale invariant instanton solutions for 
instanton size above ~ 2.3 lattice spacings. The corresponding "improved" topological charge stabilizes 
at an integer value after few cooling sweeps. At zero temperature the susceptibility calculated from this 
charge (about (195 MeV) 4 for SU(2) and (185 MeV) 4 for SU(3)) agrees very well with the phenomcno- 
logical expectation. At the minimal amount of cooling necessary to resolve the structure in terms of 
instantons and anti-instantons we observe a dense ensemble where the total number of peaks is by a 
factor 5-10 larger than the net charge. The average size observed for these peaks at zero temperature 
is about 0.4-0.45 fm for SU(2) and 0.5-0.6 fm for SU(3). The size distribution changes very little with 
further cooling, although in this process up to 90% of the peaks disappear by pair annihilation. For QCD 
we observe below T c a drastic reduction of the topological susceptibility as an effect of the dynamical 
fermions. Nevertheless also here the instantons form a dense ensemble with general characteristics similar 
to those of the quenched theory. A further drop in the susceptibility above T c is also in rough agreement 
with what has been observed for pure SU (3). We see no clear signal for dominant formation of instanton 
- anti-instanton molecules. 



1 Introduction 

The vacuum of the 4-dimensional SU(N C ) gauge theories has a non-trivial topological structure generated 
by different coverings of the spatial sphere at oo with gauge transformations from the SU(2) subgroups. 
The tunneling field configurations between these different vacua, the instantons, are self-dual solutions 
of the classical, euclidean equations of motion of finite euclidean action and integer topological charge: 
So — J d A xs(x) = 8-7T 2 , Qq = J d 4 xq(x) = ±1. The so called 't Hooft ansatz 



s(x) = F 2 (x) = - 
P 



(i) 



describes an isolated instanton of size p, centered at a; . Obviously these are scale invariant solutions. 
Superpositions of N instantons or anti-instantons also corresponds to (higher) minima of the action: 
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S = NSq, however a pair instanton - anti-instanton is not a minimum and its total action depends, 
among other, by the amount of "overlap" w = pipa/\xPj — x A \ 2 . 
Effects associated with instantons are: 

1. Ua(1) symmetry breaking. This effect is succinctly described by the Witten - Veneziano formulajlJ 
relating the topological susceptibility x and the rj' mass: 

X = ^lim ^Jdx*< T(g(x)?(0)) >= ^-(m* + m v' - 2 0- ( 2 ) 

Here x has to be obtained in a quenched simulation. Since the instantons are essentially an SU(2) 
phenomenon, the result should not depend strongly on N c . The empirical data for the RHS of (Q) give 
X = (180MeV) 4 . Notice that the presence of light fermions inhibits the topological susceptibility (the 
Dirac operator for massless fermions has zero modes for Q ^ 0). 

2. Chiral symmetry breaking. By inducing zero modes of the Dirac operator instantons may control 
the chiral symmetry breaking, due to the Banks-Casher relation < tpip >= — < p(Q) > in which p(A) is 
the spectral density of the Dirac operator. Spontaneous symmetry breaking then results from a nonzero 
density of eigenvalues at A = 0. 

3. Dynamical effects. Instantons are believed to influence the dynamics of the intermediate distances 
(at the scale of \ fm, say) and therefore the hadron properties. Their infrared effects are, however, less 
certain, in particular there is an open question whether they lead to confinement. 

The U a(X) symmetry breaking involves only the integrated correlation x- This can be tested directly in 
numerical simulations. The effects in 2. and 3. depend on details of the instanton population. Predictions 
can be made with help of models. Numerical simulations are then asked to test the ingredients on which 
these models are based. The features relevant for our discussion are: 

(i) The instantons can be in a gas, liquid, or crystalline phase depending on their density and overlap. 
The dilutcness of the ensemble is expressed by the "packing fraction" / = 7r 2 /2 < p A > N/V. 

(ii) The size distribution controls the dilutencss and the 1-R properties. In the dilute gas approximation 
the estimation of the fluctuation determinant permits to predict a power law for small sizes. In case 
of a dense population of instantons, by just integrating in a finite volume and assuming convergence of 
the thermodynamic limit one can derive some rather general features of the size distribution ||. This 
leads to a different power law behavior for the small size distribution. We can write: 

P(p)~P p , Pdiiute = -5 + b, Pdense = -f +4(6-4)/6, & = HJV c /3 (3) 

(iii) Denoting by Ni(Na) the number of instantons (anti- instantons) in a configuration, and with 
N = Ni + N A , we write: 

c^{(N*)-(Nr)/{N). (4) 

The quantity c tests the Poisson character of the N - distribution. For a poissonian distribution c = f , 
while low energy sum rules [Q suggest c = | . 

(iv) In the "standard" model || the QCD instanton ensemble below the critical temperature is an 
interacting instanton liquid with density about N/V =1 fm~ 4 . The instanton size distribution has an 
infrared cut-off and an average instanton size of ~ 1/3 fm. The dynamics of the chiral phase transition is 
driven by a rearrangement of the instanton ensemble at T c - I-A molecules are formed, with a tendency 
to orientation along the euclidean time axis (this effect is due to light fermions). 

2 Method 

Topology analyses for lattice regularized fields have the usual two sources of systematic errors: cut- 
off effects and finite size effects. Short range topological fluctuations, at the level of the cut-off, are 
unphysical. They typically have an action lower than So (dislocations) and can spoil the determination 
of the susceptibility and distort the properties of the instanton ensemble. The finite size of the lattice, on 
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Figure 1: Evolution of the size of instantons under improved cooling (left) and the (improved: 5Li or 
Wilson) action of instantons of various sizes (right). For the determination of p pea k see [6]. 



the other hand, distorts the properties of the ensemble at sizes beyond 1/2 (I: the physical lattice size). 
To cope with this latter problem we must use a large enough lattice compared with the range of distances 
we are interested in and check the dependence on the boundary conditions. The short distance problem 
is more involved. The method we have adopted for obtaining sensible densities is to cool the Monte Carlo 
configurations to smooth out the UV-noise, and in particular to eliminate the dislocations. However, the 
usual plaquette action has, strictly speaking, no instanton solutions since for them the action is always 
smaller than So and becomes lower with decreasing size. Under Wilson cooling the instantons therefore 
shrink and after some time disappear. We have used instead an improved action which has true, scale 
invariant instanton solutions for sizes p > po with a threshold pg ~ 2.3a - see Fig. |l|. This action involves 
5 Wilson loops and a similar construction yields an improved charge operator which is practically an 
integer already on rather rough configurations |(|. 

Cooling proceeds by local minimization of the action and leads to the global minimum, which is 
determined by the topological sector Q. From the point of view of our analysis it has two effects. On 
the one hand it smoothes out the short distance fluctuations, including dislocations and makes apparent 
the physical distance structure. On the other hand it annihilates I - A pairs depending on the overlap 
between the opposite charges. 

Since the (improved) charge operator Q stabilizes to an integer value within 1 - 2% already after very 
few cooling sweeps (5 - 10) the effects of type 1. above are well defined in our method without further 
tuning, monitoring, etc. The susceptibility and charge distributions show good scaling behavior which 
implies that they reveal a physical structure. The ensemble properties, on the other hand will depend, 
at least in part, on the degree of cooling. This affects the discussion of points 2. and 3. above and the 
corresponding tests for instanton models. 



3 Analysis 

We have analyzed the 5/7(2) and S/7(3) Yang-Mills theories at T = 0. We varied the cut-off and the 
lattice size such as to have the same physical volume at different discretization scales. For the SU(2) 
analysis we used twisted boundary conditions ||. We have studied QCD with two flavors of staggered 
fcrmions of mass ma = 0.008 at 3 temperatures around T c on configurations kindly provided to us by the 
M1LC collaboration. The lattices and the main results are given in Table |3|. 

TOPOLOGICAL SUSCEPTIBILITY 

As can be seen from the Fig. |^ the T = topological susceptibility is practically independent of cooling 
and scales correctly. x 1/4 shows a value of 195 - 200 MeV for 5/7(2) and 180 - 185 MeV for 5/7(3) with 
a typical error of about 5%, in very good agreement with experiment following (|J) - see also Table 0. 
This coincides with the results obtained using improved operators |8| , underrelaxed Wilson cooling 0, 
constrained smearing flo[| and with the fermionic overlap formalism flllfl (where the charge density itself 



has been tested). The results of ]12| are about 10% higher, still compatible with (g) inside 2 standard 
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Table 1: Lattices and cooling results 



Lattice 


Sw. 


(N) 


c 


(N)/fm* 


(1) SU(2),0 = 2A 
a = 0.12/m 
12 4 , 2l7conf. 


20 


10.4 


0.4 


2.4 


50 


5.0 


0.5 


1.2 


150 


2.45 


0.7 


0.6 


300 


2.12 


1.0 


0.5 


(2) SU(2),/3 = 2.6 
a = 0.06/m 
24 4 ,115con/. 


20 


63.0 


0.6 


14. 


50 


15.6 


0.65 


3.6 


150 


6.5 


0.65 


1.5 


300 


3.8 


0.5 


0.9 


(3) SU(3),/3 = 5.85 
a = 0.134/m 
12 4 , 120con/. 


20 


12.10 


0.36 


1.81 


50 


5.93 


0.5 


0.89 


150 


2.52 


1.0 


0.38 


(4) SU(S),0= 6.0 
a = 0.1/m 
16 4 , 120con/. 


20 


21.63 


0.5 


3.30 


50 


9.38 


0.5 


1.43 


100 


5.33 


0.5 


0.81 


(5) QCD,f3 = 5.65 
a = 0.115/m 
24 3 12,31con/. 


20 


50.68 


0.41 


1.75 


50 


14.90 


1. 


0.51 


100 


7.77 


0.8 


0.27 


150 


5.68 


1. 


0.20 


(6) QCD,f3 = 5.725 
a = 0.103/m 
24 3 12,43con/. 


20 


42.95 


0.57 


2.30 


50 


8.86 


0.82 


0.47 


100 


2.81 


0.98 


0.15 


150 


1.84 


1.03 


0.10 


(7) QCD,f3 = 5.85 
a = 0.0855/m 
24 3 12,20con/. 


20 


35.0 


0.4 


3.95 


50 


5.55 


0.84 


0.63 


100 


0.50 


0.9 


0.06 



deviations (see p3| for a possible explanation). We can conclude therefore that at present the T = 
topological susceptibility is under control and in good agreement with the Witten-Venetiano formula. 



Table 2: Average charge and topological susceptibility 



Lattice 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(IQI) 


1.6(1) 


1.5(2) 


1.8(2) 


1.7(2) 


2.0(4) 


0.9(1) 





X 1/4 MeV 


195(4) 


194(14) 


184(6) 


182(7) 


134(10) 


102(5) 






The finite temperature deconfining transition of the pure Yang-Mills theory has been shown to pro- 
duce a strong drop in the susceptibility j|, [To). We have analyzed here the finite temperature QCD. Below 
T c we found as expected a reduction in the susceptibility (by a factor 3), indicative of the dynamical 
fcrmion effects. Above T c (which is at ~ WQMeV for the MILC simulation) we found a pronounced 
drop in x, similar to the quenched case. As already observed before for other staggered fermion simu- 
lations jlj] also the MILC configurations show long range correlations in Q and explore few topological 
sectors. The results quoted in Table [| are obtained after symmetrizing per hand the charge distribution 
(since this is a symmetry of the action), the errors given are only the statistical ones. We have noticed 
such metastabilities also in the pure SU(3) simulation, there, however, one can afford to wait until the 
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Figure 2: Topological susceptibility for quenched SU{2) and SU(3). 




Figure 3: Isosurfaces for the charge density of a typical configuration of lattice (6), Table 1., after 20 
sweeps of improved cooling; the two circles indicate a "molecule" (left) and a "caloron" (right). 



configurations decorrelate. 

THE INSTANTON ENSEMBLE 

To observe not merely the global charge of a configuration but its structure describable in terms of instan- 
tons and anti-instantons a certain degree of smoothing must be attained (see Fig. || for an illustration). 
Using our improved cooling we achieve this smoothness with a number of cooling sweeps between 20 and 
50, depending on the cut-off a. We did not try to tune a rescaling law for the amount of cooling 0, since 
for many of the results this was not necessary due to inherent stability properties of our cooling. 

As can be seen from the tables, when the configurations are just smooth enough for us to see instantons 
and anti-instantons we typically have a ratio N/ \Q\ ~ 6 — 10. This ratio is reduced to Nf \Q\ ~ 1 — 2 in the 
further cooling by pair annihilation, leaving only the (anti-)instantons corresponding to the topological 
sector. The first remark is that the size distribution shows only little variation during this process - see 
Fig. ^ for the quenched simulation. This implies that all the (anti-)instantons which we observe from 
the start obey the same distribution. The average sizes are in the region of 0.4fm for SU(2) and 0.6fm 
for SU(3) and QCD. The small size part of the distribution departs from the dilute gas power law and 
shows tendency toward agreement with ||, [l^] (p ~ 0.8 for SU{2) and 1.5 for SZ7(3) - compare eq. (0)). 

The density of (anti-)instantons is high to start with. Because of this high density and of the rather 
large sizes, the packing fraction here is also high. The overall distribution is here non-poissonian and in 
agreement with the low energy sum rules, see eq. (^) and Table || (while it approaches the poissonian 
distribution for instantons or anti-instantons after long cooling) . This suggests a dense ensemble with a 
spatially random distribution biased toward small \Q\. Similar results at T = have been reported in 

§■ 

To the extent we could follow the structure described by this ensemble back to shorter cooling, it does 
not change in character but becomes more and more obscured by the accumulation of short range ripples 
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and dislocations. This corroborates with the stability of distributions under further cooling (Fig. |j) and 
their correct scaling behavior to suggest the physical relevance of this ensemble. Since the short cooling 
needed before we can observe the instantons and anti-instantons might have also annihilated some pairs 
in the physically relevant size region, the figures we obtain for the density of (anti-)instantons are lower 
bounds. 

For QCD at high T molecule formation with preferred orientation should be revealed by an anisotropy 
in the density - density correlations along space and along (euclidean) time directions. The anisotropy we 
observe at very short cooling is compatible with zero at the level of two standard deviations and very low 
in relative value (compared to the correlation itself) This signal is too weak to support the scenario 
of Q, however the quark mass might not be yet small enough to promote a strong effect. See also Fig. |[ 

We conclude that we observe a dense ensemble randomly distributed in space, but with a bias toward 
small \Q\. This ensemble is revealed after a short cooling which essentially has filtered out dislocations 
and strongly overlapping I A pairs. We remark as a question of principle that any separation 

q(x) — instantons + trivial fluctuations 
is inherently ambiguous in the case of such strongly overlapping pairs. It is therefore not clear to us 
whether this latter structure is reasonably described as instantons and anti-instantons or a description 
in terms of other kind of excitations is more suitable for it. 
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